Magnetic-Field Dependence of the YbRh 2 Si2 Fermi Surface 
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"Large" Fermi surface: YbRh 2 Si 2 LDA+SOC 
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Magnetic-field-induced changes of the Fermi surface play a central role in theories of the exotic 
quantum criticality of YbRli2Si2. We have carried out de Haas-van Alphen measurements in the 
magnetic-field range 8 T < H < 16 T, and directly observe field dependence of the extremal Fermi 
surface areas. Our data support the theory that a low-field "large" Fermi surface, including the Yb 
4/ quasihole, is increasingly spin split until a majority-spin branch undergoes a Lifshitz transition 
and disappears at .Ho ~ 10 T, without requiring 4/ localization at Hq. 

PACS numbers: 71.18.+y, 71.20.Eh, 71.27. +a, 75.20.Hr 



The heavy fermion compound YbRIi2Si2 has been the 
focus of much recent interest because it can be tuned 
through a field-induced quantum critical point (QCP) 
which is both easily accessible to experiment and ap- 
pears to exhibit a new class of "local" quantum criti- 
cality [ij. It has been argued that this scenario involves 
dramatic Fermi surface (FS) changes as a function of tem- 
perature and magnetic field [l|, [2J, |3|. In this Letter, we 
have used de Haas-van Alphen (dHvA) measurements to 
probe these proposed changes. 

The main focus of attention in YbRli2Si2 has so far 
been at low magnetic fields, where the magnitude of the 
Hall coefficient shows a series of crossovers, extrapolating 
to a discontinuous jump at the H c \\c « 0.6 T quantum 
critical point. This has been interpreted as signalling a 
sudden Fermi surface reconstruction from a small to a 
large FS [4| ; such an interpretation is controversial, how- 
ever, due to the sensitivity of the YbRli2Si2 Hall effect 
to small changes in /-occupancy [a, [6j and/or changes in 
quasiparticle scattering as antiferromagnetic fluctuations 
give way to ferromagnetic fluctuations for H > H c [7| . 

At higher fields, specifically at HqJ-c w 10 T, 
crossovers of dc-magnetization, T 2 coefficient of resistiv- 
ity, specific heat, and linear magnetostriction coefficient 
have been interpreted as evidence for a second Fermi sur- 
face reconstruction, back to the small FS, as the Zeeman 
energy becomes comparable to the Kondo energy scale. 
This has been suggested to be a continuous change Q or 
a transition [2J of the total Fermi volume. 

The primary goal of our work was to look near H 
for this second proposed Fermi surface change. Because 
dHvA is limited to high fields, we could not directly ac- 
cess the low-field QCP, but the Fermi surface behavior 
near Hq sheds light on the electronic structure between 
H c and Hq, and therefore provides experimental con- 
straints for theories of YbRli2Si2 quantum criticality. 

The central issue with regard to the proposed Fermi 
surface reconstructions discussed above is whether or not 
the Yb 4/ quasihole is included in the Fermi volume. 
In the "small" Fermi surface case (also known as "4/ 
localized" or "Yb 3+ "), it does not contribute, whereas 
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FIG. 1: (color online) "Small" and "large" Fermi surfaces 
calculated with the LDA + spin-orbit coupling method. Fol- 
lowing convention from Ref. [ill ], the sheets are labeled donut 
U D" jungle-gym "J," and pillbox "P." From an electron point 
of view, the dark (blue online) side of each sheet is the occu- 
pied side, and the light (yellow online) side is the unoccupied 
side, such that D is a hole surface and P is an electron sur- 
face. Note that D is a torus in the small Fermi surface case 
only. 



in the "large" Fermi surface case (also known as "4/ 
itinerant"), it does. To illustrate these two cases, we 
have performed all-itinerant LDA + spin-orbit coupling 
calculations § for LuRh 2 Si 2 (small FS) and YbRh 2 Si 2 
(large FS) using the wien2k density functional theory 
code [9j . The large-scale features of our calculated Fermi 
surfaces, visualized in Fig. [1] via the xcrysden pro- 
gram [id ], are in good agreement with those published 
previously @, EH, E3, Lll • 

In order to probe the Fermi surface of YbRli2Si2 as 
the field is increased through Hq, we have conducted 
ambient-pressure dHvA measurements from 8-16 T and 
30-600 mK, using the standard field modulation method. 
The high quality single-crystal samples studied were 
~0.2 x 1 x 2 mm platelets grown from indium flux [l3j |. 
with a residual resistivity ratio of about 100. The mag- 
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FIG. 2: (color online) YbRli2Si2 quantum oscillations in the 
14-16 T range, with field aligned ~7° from (110) toward 
(001). A Fast Fourier Transform (FFT) in 1/H of the dis- 
played data is inset. 



netic field was rotated across a range of 45° from a = 
(100) to (110) and 60° from (110) toward c = (001). Fig- 
ure [2] shows a sample trace of the as-measured dHvA 
oscillations and an associated Fast Fourier Transform 
(FFT) in the 14-16 T field range. 

Figure [3] shows the field- angle dependence of our mea- 
sured quantum oscillation frequencies in the 14-16 T field 
range, below 50 mK. For comparison, previously pub- 
lished data from torque measurements on samples grown 
in the same facility [13J|, as well as predicted dHvA fre- 
quencies from the calculated small and large Fermi sur- 
faces of Fig. [TJ are included. Predicted dHvA frequencies 
and band masses were extracted from each calculated 
Fermi surface sheet using an algorithm described else- 
where [l4|. With the exception of the extra frequency 
near 6 kT beyond -15° in the (HO)-(OOl) plane, and 



persistence of observed frequencies over a larger angular 
range in this plane, there is reasonable qualitative agree- 
ment between the small Fermi surface calculation and 
experiment. Quantitative agreement is generally not ob- 
tained for heavy fermion compounds since the flatness of 
the bands makes the FS size depend sensitively on small 
shifts of the Fermi energy. The qualitative agreement 
between the large Fermi surface calculation and experi- 
ment is significantly poorer, since the calculation misses 
numerous branches in the (100)-(110) plane. 

To highlight the correspondence between experiment 
and theory, matching pairs of orbits have been labeled in 
Fig. [3] with the same letter and number combination — 
uppercase letters for calculated orbits from the small FS 
case and lowercase letters for experimental data. Pre- 
dicted major orbits on the small J surface are at higher 
or lower frequencies and are not shown in Fig. [3J as we 
were unable to observe them: all of our experimentally 
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FIG. 3: (color online) dHvA frequency vs magnetic-field an- 
gle. Closed (blue online) circles depict the current experimen- 
tal data (14-16 T field range), open squares show previously- 
published data (12-28 T field range) [igj|. and lines are from 
the calculated small (solid, green online) and large (dashed, 
red online) Fermi surfaces shown in Fig. [l] Calculated or- 
bits from the small FS case are labeled with uppercase (green 
online) letters and numbers, whereas corresponding experi- 
mental orbits are labeled with lowercase (blue online) letters 
and numbers. 



TABLE I: Calculated band masses (mb) from the small Fermi 
surface in Fig. [T] measured effective masses (to*) from the 
dHvA experiment (14-16 T field range), and resulting mass 
enhancements (m* /nib), for calculated and measured orbits 
in Fig. El 



LDA + SOC orbit 


Expt. orbit 


m b (m e ) 


m* (m e ) 


m* /m,b 


Dl 


dl 


1.00(6) 


10.1(5) 


10.1(8) 


D-l 
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0.82(2) 


8.6(5) 


10.5(7) 


D5 


db 


0.71(3) 


8.6(5) 


12.1(9) 
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0.62(5) 


6.8(3) 


11(1) 


D7 
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0.72(1) 


10(1) 


14(1) 


D8 


(18 


0.554(4) 


5.00(5) 


9.0(1) 


Dll 


dll 


1.69(9) 


9.02(9) 


5.3(3) 



detected frequencies correspond to hole orbits on the D 
sheet. Calculated band masses, measured cyclotron ef- 
fective masses, and corresponding mass enhancements for 
measured orbits and their theoretical counterparts are 
listed in TablelU the average mass enhancement is 10(1). 
Within the 0.1 kT FFT resolution of our experiment, no 
spin splitting of the dHvA frequencies is observed, and 
the temperature dependence of all measured oscillations 
is well described by the standard Lifshitz-Kosevich rela- 
tion. 

As our aim was to investigate possible Fermi surface 
changes across Hq, it is important that we were able to 
follow the strongest oscillations (d5, d6, d8 and dll from 
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FIG. 4: (color online) Measured dHvA frequencies as a func- 
tion of magnetic-field strength, normalized to each orbit's fre- 
quency at 15.5 T. Orbit labels correspond to those in Fig. [3] 
the high-field frequencies for each orbit are: Fds — 5.4 kT, 
F dG = 3.5 kT, F ds = 2.6 kT, and F dn = 6.1 kT. The inset 
schematically shows the correspondence between a true Fermi 
volume change (i) and that observed via back-projection in 
a dHvA experiment (ii). The frequencies we follow lowest in 
field show the S-shape of inset (ii) . 



Fig. El) from 16 T (well above H ) to 8 T (well below 
Hq). Figure [4] shows the measured frequency of these 
oscillations as a function of field strength, normalized 
to the value at 15.5 T. Remarkably, all of the four D- 
surface orbits that we can follow below Hq exhibit the 
same behavior: roughly constant above about 11 T, rising 
sharply as the field is reduced through Hq, and perhaps 
tending to saturate below ~9 T. 

Interpretation of the field dependence of a measured 
dHvA frequency, f m , requires some care, because f m is 
not directly related to the "true" frequency, ft , but rather 
it is the back-projection to H = of the line tangent to 
ft, as illustrated in the inset of Fig. |4] [15]. For example, 
at a particular field H A , f mjA = f t>A - H A (df t /dH) A . 
Observed field dependence of f m usually reflects changes 
oidft/dH rather than of ft itself. In YbRli2Si2, we know 
that the magnetization roughly saturates above Hq [2| , so 
the slope \dft/dH\ should be lower for H > Hq. Thus, 
the behavior must be roughly as shown in the inset of 
Fig. 21 with a gradual change in f m corresponding to a 
gradual change in f t and therefore the Fermi surface. We 
conclude that we are observing a continuous shrinking of 
the measured £)-surface at a rate that slows above Hq. 

It is tempting to infer that the shrinking D surface 
means that the Yb 4/ quasihole is gradually disappear- 
ing from the Fermi volume.and that Hq represents the 
completion of this process [3(. However, the traditional 
Luttinger's theorem-based view is that at T = K con- 
tinuous changes in total Fermi volume should not occur. 



If the Yb 4/ quasihole were to localize across Hq , a dis- 
continuous jump in dHvA frequencies would be observed, 
as in CeRhIn 5 under pressure [16J . Since the observed fre- 
quencies change continuously across Hq, 4/ localization 
must not occur at this field. 

An alternative scenario for the behavior at Hq has been 
developed by Kusminskiy et al. [17J]. They consider the 
spin splitting of the Fermi surface at high fields, and 
propose that Hq represents a Lifshitz transition, where 
a heavy, majority-spin branch of the large Fermi sur- 
face disappears, leaving only a moderate-effective-mass 
minority-spin branch of the large Fermi surface. The 
4/ localization transition is predicted to occur at much 
higher fields. This phenomenon occurs in both static [17j 
and dynamic [l8j mean field theory treatments. 

The Kusminskiy scenario was developed for a toy, one- 
band model, but we show in Fig. O how it might work for 
the actual band structure of YbRli2Si2, using a simplified 
sketch of the D and J energy bands along the Z-U direc- 
tion in the Brillouin zone. The bands giving rise to the 
small Fermi surface [solid grey lines in Fig.[5ja)] hybridize 
with a virtual /-level near the Fermi energy, Ep. The 
resulting many-body quasiparticle bands [dashed lines in 
Fig. [DJa)] have a large Fermi surface (i.e., the /-hole is 
now included in the Fermi volume), normally assumed to 
be the same as the large LDA Fermi surface [dotted grey 
lines in Fig. E{a)]. When the magnetic field polarizes 
the quasiparticle bands [dashed lines in Fig. Ob)], the 
majority- and minority-spin branches split: the minority 
hole quasiparticle _D-band sinks, such that its Fermi sur- 
face resembles the small D FS sheet; the majority hole 
quasiparticle £>-band rises to resemble the small J FS 
sheet, where it becomes the opposite-spin counterpart of 
the minority hole quasiparticle J-band; finally, the ma- 
jority hole quasiparticle J-band rises, and Hq is the field 
at which the quasiparticle Fermi surface corresponding 
to this band has grown to encompass the entire Brillouin 
zone and disappears. Thus, through this process, it is 
possible for the minority branch of the large D Fermi 
surface to resemble the small D Fermi surface without 
requiring 4/ localization at Hq. 

The situation is reminiscent of the heavy fermion com- 
pound CeRu2Si2, which exhibits a metamagnetic tran- 
sition (MMT) at applied fields of H = (7.7/ cos H ) T, 
where Oh is the angle between the field direction and the 
crystallographic c-axis. dHvA measurements below and 
above the MMT were matched to band structure predic- 
tions for large and small Fermi surface cases respectively, 
implying a sudden localization of the Ce 4/ state at the 
MMT [19(. However, as in YbRli2Si2, measurements of 
other CeRu2Si2 bulk properties rule out the possibility 
of a first order transition at the MMT, leading to the 
proposal of a continuous Fermi surface Lifshitz transi- 
201. 
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tion 

The modified Kusminskiy model discussed above is in 
accord with both Luttinger's theorem and our experi- 
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FIG. 5: (color online) A schematic representation of the 
YbRli2Si2 bands near the Fermi energy, Ef, following Kus- 
minskiy et al. [17| . Bands associated with the "D" and "J" 
Fermi surface sheets of Fig.fJJare labeled. Panel (a) shows the 
band structure without spin splitting: solid gray lines for the 
small Fermi surface case, and dashed (blue, red and green on- 
line) lines for the many-body quasiparticle bands whose Fermi 
surface coincides with the "large" LDA FS (dotted lines). The 
virtual /-hole is shown by a gray horizontal line in panel (c). 
Panel (b) adds spin splitting to the quasiparticle bands, caus- 
ing one minority-spin surface to resemble the unsplit small D 
surface, and the largest majority-spin band to no longer cross 
Ef. Panels (c) and (d) show magnified views of panels (a) 
and (b), respectively, near Ef- 



mental results. A key point is that we do not resolve 
spin splitting of the orbits in any of our measurements, 
suggesting that what we observe are indeed small-D-like 
oscillations that arise from one spin direction only, as in 
Fig. [3b). If the 4/ quasihole had localized at Ho, all 
of the observed orbits would be spin split. In contrast, 
in the Kusminskiy model, the majority-spin quasiparti- 
cle branch has a completely different topology from the 
minority-spin branch. As the minority-spin quasiparti- 
cle branch shrinks to resemble the small D LDA Fermi 
surface, it becomes a topological torus, allowing the ap- 
pearance of orbits D6, 7, 8, 9, 10, and 12, which thread 
through the center of the torus. The corresponding 
majority-spin branch is never toroidal, and thus would 
not have these orbits. Moreover, as this majority surface 



grows to resemble the small J LDA Fermi surface, its 
topology will radically change, such that orbits Dl, 2, 3, 
4, and 11 also vanish. 

We therefore believe that our result, combined with 
Luttinger's theorem, rules out the 4/ localization sce- 
nario and supports the approach of Kusminskiy et 
al. [17J : the high-field Fermi surface is a spin split version 
of the large YbRli2Si2 Fermi surface, rather than a direct 
realization of the small Fermi surface. This, in turn, im- 
plies that the Fermi surface includes the 4/ degrees of 
freedom for H c < H < Hq. 
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